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1.  Introduction 


While  most  projectiles  in  use  today  behave  as  rigid  bodies  while  in  flight,  a  notable  number  of 
projectiles  are  purposefully  designed  to  carry  a  liquid  payload.  For  example,  smoke  screens 
delivered  by  artillery  rounds  consist  of  a  typical  spin-stabilized  shell  containing  a  canister  filled 
with  white  phosphorous  {1-4).  Another  example  is  new  less-than-lethal  projectiles  having  a 
concentric  cylindrical  cavity  filled  with  liquid  that  delivers  this  payload  to  a  target  upon  impact 
(5).  Finally,  some  projectiles  are  designed  to  be  general  payload  delivery  shells,  including 
delivery  of  medical  supplies,  such  as  intravenous  fluid  bags  (6). 

There  can  be  a  significant  difference  in  flight  behavior  between  liquid-  and  solid-filled 
projectiles.  The  difference  is  caused  by  motion  of  the  liquid  inside  the  spinning  projectile.  This 
motion  causes  forces  to  act  on  the  projectile  body,  which  can  prematurely  terminate  the  flight  by 
instability.  Characteristics  of  this  instability  are  sharp  increases  in  angle  of  attack  accompanied 
by  large  changes  in  spin  rate  {7-10).  At  gun  launch,  the  motion  of  the  projectile  causes  the  fluid 
to  spin  up  in  a  time-dependent  manner,  but  it  subsequently  achieves  steady  flow. 

Prediction  of  instability  induced  by  a  liquid  payload  installed  in  a  coning  projectile  has  been 
analyzed  by  computational  fluid  dynamic  (CFD)  simulations  and  with  analytic  linear  fluid 
theories  based  on  spatial-eigenvalue  methods.  The  CFD  solutions  are  generally  best  suited  to 
liquids  with  a  low  Reynolds  number  (Re),  while  spatial-eigenvalue  methods  can  encompass  a 
much  broader  range  of  Re.  A  good  survey  of  liquid-filled  projectiles  with  a  focus  on  flight 
instabilities  is  given  in  Sedney  {11).  Generally,  analytic  theories  are  a  composition  of  inviscid 
flow  and  viscous  boundary-layer  corrections  {12-14).  The  prediction  of  flight  stability  of  a 
liquid-filled  projectile  has  also  been  studied  using  tricyclic  linear-projectile  theory  {15).  This 
analysis  assumes  the  effect  of  a  liquid  payload  is  similar  to  the  Magnus  effect.  Spectral  analysis 
used  to  numerically  compute  liquid-fill  induced  moments  has  been  incorporated  into  six-degree- 
of- freedom  (6-DOF)  simulations  {16).  The  angular  motion  of  a  low-Re  liquid-filled  projectile 
has  also  been  simulated  using  a  precomputed  table  of  liquid  moments,  obtained  from  CFD,  in 
6-DOF  calculations  {1 7). 

The  present  report  meshes  a  well-developed  spatial-eigenvalue  theory  directly  into  a  standard 
6-DOF  projectile  flight  dynamic  model.  Hence,  a  well-established  rotating  liquid  model 
calculates  liquid  moments  at  each  time  step  of  the  6-DOF  trajectory  numerical  simulation 
without  the  need  of  interpolation  to  gain  the  effects  of  liquid  payloads.  This  report  further 
bridges  the  gap  between  the  body  of  literature  on  effects  of  viscous  liquid  payloads  on  projectiles 
and  6-DOF  projectile  flight  dynamic  modeling.  The  report  begins  with  a  review  of  rigid 
projectile  flight  dynamic  modeling  along  with  a  description  of  modeling  a  rotating  liquid  in  a 
cylindrical  cavity.  The  two  models  are  subsequently  integrated  such  that  a  projectile  flight 
dynamic  model  with  a  liquid  payload  results.  To  highlight  the  utility  of  the  methodology,  the 
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flight  dynamic  model  is  exercised  on  an  example  shell.  Comparisons  are  made  between  a  liquid- 
filled  projectile  and  a  similar  solid-filled  projectile  highlighting  the  predictive  capability  of  the 
new  model. 


2.  Projectile  Flight  Dynamic  Model  With  a  Liquid  Payload 


A  typical  6-DOF  rigid  projectile  model  is  employed  to  predict  the  dynamics  of  a  projectile  in 
flight.  These  equations  assume  a  flat  Earth.  The  well-known  6-DOF  states  comprise  the  three 
translational  components  describing  the  position  of  the  projectile’s  center  of  mass  and  the  three 
Euler  angles  describing  the  orientation  of  the  projectile  with  respect  to  the  Earth.  Figures  1  and  2 
provide  a  visualization  of  the  degrees  of  freedom. 


Figure  1.  Projectile  position  coordinate  definitions. 
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Figure  2.  Projectile  orientation  definitions. 


The  equations  of  motion  (18)  derived  in  the  no-roll  frame  are  shown  in  equations  1-4. 
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Note  that  the  following  shorthand  notation  for  trigonometric  functions  is  used  in  equations  1—4: 
=  sin(a),  ca  =  cos(a),  and  ta  =  tan(a). 
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The  force  acting  on  the  projectile  in  equation  3  comprises  the  weight  force  (W)  and  the 
aerodynamic  force.  The  aerodynamic  force  is  split  into  a  standard  force  (A)  and  a  Magnus  force 
(M).  The  combination  of  forces  is  expressed  in  equation  5. 
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Equation  6  provides  the  expression  for  the  weight  force  in  the  no-roll  coordinate  system. 
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Equation  7  provides  the  expression  for  the  aerodynamic  force  in  the  no-roll  coordinate  system. 
This  force  acts  upon  the  projectile  at  the  aerodynamic  center  of  pressure. 
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Equation  8  provides  the  expression  for  the  Magnus  force  in  the  no-roll  coordinate  system.  The 
Magnus  force  acts  upon  the  projectile  at  the  Magnus  force  center  of  pressure. 
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Equations  7  and  8  are  based  on  Mach  number-dependent  coefficients  and  the  total  aerodynamic 
velocity  given  in  equation  9. 

V  =  Vu"  +  v“  +  w"  .  (9) 

The  moment  acting  on  the  projectile  in  equation  4  comprises  the  moment  due  to  the  standard 
aerodynamic  force  (A),  the  moment  due  to  the  Magnus  aerodynamic  force  (M),  the  unsteady 
aerodynamic  moment  (EIA),  and  the  liquid  payload  moment  (L)  as  shown  in  equation  10. 
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The  moment  due  to  the  aerodynamic  force  is  expressed  in  equation  1 1 . 
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The  moment  due  to  the  Magnus  force  is  expressed  in  equation  12. 
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The  unsteady  aerodynamic  moments  acting  on  the  projectile  are  expressed  in  equation  13. 
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The  coefficients  used  in  this  model  are  specific  functions  of  the  Mach  number  of  the  projectile. 
For  fin-stabilized  projectiles,  Magnus  force  and  moment  are  typically  ignored  since  its  effect  is 
rather  small  for  slowly  rolling  projectiles. 

The  angular  motion  of  a  projectile  is  altered  by  the  inertial  waves  propagating  through  the  liquid 
payload.  In  turn,  these  waves  act  on  the  walls  of  the  liquid  cavity  and  generate  a  liquid  moment. 
The  impact  of  this  moment  can  have  a  devastating  impact  on  the  projectile’s  angular  motion 
(19-21).  For  this  report,  the  liquid  moment  is  calculated  from  steady  linear-projectile  theory, 
which  incorporates  solutions  of  the  linearized  Navier-Stokes  equations  (18).  In  general,  this 
liquid  moment  model  is  a  linear  combination  of  fast  and  slow  mode  liquid  moment  contributions. 
In  section  2  of  his  report,  Murphy  (19)  states  that  for  a  steady-state  coning  projectile  with  a 
liquid  payload,  .  .  only  the  fast  mode  motion  is  adversely  affected  by  the  liquid  side  moment.” 
For  this  reason,  only  the  fast  coning  mode  is  used  to  calculate  the  liquid  moment  in  this  paper. 
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Figure  3  shows  the  orientation  of  body  unit  vectors  (lB  ,  JB  ,  KB  )  to  velocity  unit  vectors 
(lv  ,  Jv  ,  Kv  ) .  Appling  the  liquid  theory  (18)  to  figure  3  allows  the  liquid  transverse  and  roll 
liquid  moments  to  be  written  as  shown  in  equation  14. 


Figure  3.  Liquid  transverse  velocity  orientation 
definitions. 
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(14) 


The  moment  coefficients  Clsm  and  Clim  depend  on  the  fast  mode  coning  frequency  T,  fast  mode 
undamping  s,  the  cavity  aspect  ratio  c/a,  liquid  Re,  and  the  magnitude  of  axial  spin  (19).  Values 
of  the  fast  mode  pair  (T,  s)  are  obtained  from  equation  15. 


T  = 


2p2  -71PCnad2IyRcxV2  |  Ix+a2mLCL1M 


2x/2pIY 

D(CMQT+CYPARMX)  Cna(IyT-Ix) 
2  M 


TipD2 


V 


8  = 


8pT(2IYT-Ix) 


-  + 


2L 


a  mLCLSM 
(2IyT-Ix) 


(15) 


6 


The  liquid  moment  coefficients  Clsm  and  Clim  are  calculated  using  the  techniques  given  by 
Murphy  (19). 


3.  Results 


The  example  simulations  given  here  consider  liquid  payloads  in  a  typical  projectile  with  launch 
properties  given  in  table  1 . 


Table  1.  Launch  properties  of  projectile. 


Physical  Parameters 

Magnitude 

lx 

0.1157  slug  ft2 

Iy,  Iz 

1.1972  slug  ft2 

m 

103.0  lb/ft 

D 

0.510  ft 

Launch  velocity,  V 

2460.0  ft  s  1 

Launch  axial  spin  rate,  p 

1500.0  s-1 

The  liquid  density  pLis  taken  to  be  3.5  times  the  density  of  liquid  water,  and  the  viscosity  n  is 
selected  such  that  launch  Re  >  8  x  106.  This  is  sufficient  to  ensure  that  the  Re  remains  large  so 
boundary  layer  analysis  (18)  adequately  governs  the  liquid  physics  for  the  entire  trajectory.  The 
liquid  cavity  is  a  cylinder  with  aspect  ratio  c/a  =  3.75  completely  filled  with  this  low-viscosity 
liquid.  Figure  4  shows  the  range  of  nondimensional  coning  frequencies,  T,  for  a  typical 
trajectory  of  this  projectile  housing  a  frozen  (solid)  liquid  payload. 
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The  Re  for  a  given  aspect  ratio  can  cause  large  variations  in  the  liquid  moment  when  subjected  to 
a  range  of  nondimensional  coning  frequencies.  Applying  steady-state  linear  liquid  theory  to  a 
payload  configuration  undergoing  coning  motion  reveals  in  figure  5  the  side  moment  coefficient 
Clsm  with  c/a  =  3.75  for  two  numbers  Re  =  8  x  1  ()':  and  8  x  10s .  These  results  depict  important 
CLSM  behavior,  where  the  peaks  indicate  a  potential  problem  due  to  large  liquid  moments  when 
the  nondimensional  coning  rate  T  «  0.088  .  Note  the  amplitude  of  the  peak  has  a  strong 
dependence  on  Re.  Peak  values  of  Clsm  are  nearly  2  orders  of  magnitude  larger  than  the 
nonpeak  values  of  Clsm- 


Figure  5.  CLsm  vs.  coning  frequency  and  Re,  resonance  T  =  0.088. 

Figure  5  suggests  that  for  a  sufficiently  large  Re,  a  free-flight  missile  having  nondimensional 
coning  frequencies  in  the  neighborhood  of  T  =  0.088  will  experience  a  significantly  larger  liquid 
moment  than  it  does  for  frequencies  outside  this  neighborhood.  Such  results  are  characteristic  of 
the  root  cause  of  projectile  instabilities  due  to  liquid  payloads.  The  aspect  ratio,  coupled  with  a 
large  enough  Re,  forces  inertial  waves  in  the  coning  fluid  to  generate  large  CLsm  coefficients 
(19).  In  general,  increasing  Re  causes  the  liquid  side  moment  to  increase  near  T  =  0.088,  and 
this  pronounced  peak  signifies  a  possible  resonant  frequency.  This  is  one  indicator  of  a  liquid 
payload  possibly  causing  the  fast  coning  frequencies  to  change  such  that  an  instability  occurs 
near  T  =  0.088.  Note  the  range  of  frequencies  given  in  figure  5  is  close  to  the  fast  mode  T 
frequencies  for  the  frozen  liquid  found  in  figure  4. 

Next  a  series  of  results  showing  the  effect  a  liquid  payload  has  on  projectile  angular  motion 
during  flight  is  shown.  Results  presented  here  are  for  a  launch  Re  =  8  x  1 06,  which  is  large 
enough  so  linear  liquid  theory  with  boundary  layer  corrections  is  valid  throughout  the  entire 
trajectory.  Figures  6  and  7  compare  roll  moment  and  spin  rate  results  for  a  frozen  liquid  and  a 
flowing  liquid.  These  variables  are  selected  because  the  rotational  physics  of  a  projectile  causing 
a  rapid  decrease  in  spin  rate  is  a  strong  indicator  of  flight  instability  caused  by  liquid  payloads 
(18). 
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Figure  6.  Projectile  roll  moment  for  frozen  and  flowing  liquid  payloads. 


The  liquid  payloads  for  the  launch  Re  chosen  in  these  examples  has  small  effects  on  liquid  roll 
moment,  as  shown  in  figure  6.  As  expected,  figure  7  shows  the  projectile  roll  rates  also  change 
slightly  compared  to  the  frozen  liquid  payload. 


Figure  7.  Projectile  spin  rate  for  frozen  and  flowing  liquid  payloads. 
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Evidently  the  increased  magnitude  of  the  projectile  roll  moment  due  to  the  flowing  liquid  is  not 
sufficient  to  cause  flight  instability  since  the  roll  performance  of  the  projectile  changes  very  little. 

The  liquid  side  moment  coefficient,  Clsm,  is  an  important  parameter  when  evaluating  the  impact 
a  liquid  payload  has  on  projectile  angular  motion.  Figure  8  shows  the  fast  mode  liquid  moment 
coefficient  as  a  function  of  T.  This  chart  shows  that  Clsm  changes  rapidly  for  the  region 
0.086  <  T  <  0.090 ,  and  the  magnitudes  of  Clsm  exhibits  the  approximation  CLSM  ~  VRe  (18) 
when  the  coning  rate  is  in  the  neighborhood  of  T  =  0.088. 


Figure  8.  Liquid  moment  coefficients  vs.  nondimensional  fast  mode  frequency. 

Figure  9  compares  the  angle  of  attack  (AO A)  for  identical  projectiles  with  payloads  supporting 
inertial  waves  and  frozen  liquid  payloads. 


Figure  9.  Comparison  of  AOA  of  frozen  and  liquid  payloads. 
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Apparently  the  inertial  wave  motion  slightly  increases  the  AOA  for  a  liquid  with  launch  Re  of 
8.0  x  106  or  8.0  x  107  compared  to  the  frozen  liquid,  but  this  increase  is  not  enough  to  cause 
concern  of  flight  instability.  In  fact,  the  results  presented  so  far  show  no  indication  of  flight 
instability  even  when  the  projectile  has  a  coning  rate  near  T  =  0.088  for  the  potential  problem  in 
Clsm  found  in  figures  5  and  8. 

In  order  to  understand  the  dynamics  of  a  liquid-filled  projectile  exhibiting  angular  instability,  we 
consider  figure  4  and  exploit  the  peak  value  of  Clsm  near  T  =  0.088.  Thus,  to  substantially 
increase  the  size  of  the  liquid  moment,  we  assume  the  payload  comprises  a  hypothetical  liquid 
with  a  decreased  viscosity  so  that  the  launch  Re  is  large  enough  to  cause  substantial  increases  in 
liquid  side  moment  during  flight.  In  particular,  we  want  this  liquid  to  generate  large  enough 
values  of  Clsm  for  T  ~  0.088  to  alter  the  angular  motion  of  the  projectile.  A  sufficiently  large 
decrease  in  projectile  spin  rate  along  with  a  large  increase  in  AOA  suggests  the  hypothetical 
liquid  may  cause  premature  termination  of  flight.  Figure  10  shows  the  fast  nondimensional  T 
dependence  of  the  side  moment  coefficient,  Clsm,  rapidly  increasing  near  T  =  0.088  for  launch 
Re  =  4  x  10  .  This  indicates  a  possible  increase  in  the  magnitude  of  projectile  angular  motion, 
which  is  consistent  with  the  resonance  configuration  presented  in  figure  5. 


Figure  10.  Comparison  of  projectile  roll  moment  vs.  payload  launch  Re. 


Figure  1 1  compares  the  AOA  of  the  projectile  with  liquid  payloads  having  three  different  launch 
Re’s. 
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The  projectile  with  launch  Re  =  4  x  10  causes  the  fast  frequency  at  t  ~  4.5  s,  T  ~  0.088,  to 
generate  an  AO  A  >50°  for  which  the  judgment  was  made  to  terminate  the  numerical  integration. 
Figure  12  shows  the  corresponding  projectile  roll  moment  L  rapidly  increasing  in  magnitude  due 
to  the  larger  Re.  Note  that  in  the  neighborhood  of  T  =  0.088,  the  fast  mode  liquid  motion  is 
responsible  for  increasing  the  magnitude  of  L. 


Projectile  Roll  Moment 


Figure  12.  Comparison  of  projectile  roll  moment  vs.  launch  Re. 
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Figure  13  shows  the  corresponding  projectile  roll  rates  for  the  launch  Re  =  8  x  10  and  4x10, 
which  correlates  with  the  size  of  liquid  moment  coefficient  CLrm  as  shown  in  equation  15.  The 

•*ll  8 

increase  in  L  for  Re  =  4  x  10  causes  the  rapid  decrease  in  roll  rate. 


o 

Figure  14  compares  corresponding  pitch  rates  of  the  frozen  liquid  and  launch  Re  =  4  x  10‘ 
configurations.  The  rapid  increase  in  q  values  near  t  ~  4.5  s  is  caused  by  the  fast  mode  liquid 
side  moment. 
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A  contrast  in  projectile  yaw  rates  for  frozen  and  liquid  payload  is  given  in  figure  15.  Again  the 
rapid  increase  in  yaw  rate  takes  place  near  t  ~  4.5  s  in  the  neighborhood  of  the  fast  mode 
resonance. 


These  calculations  show  a  significant  difference  in  projectile  angular  motion  when  the  liquid 

8 

payload  has  a  launch  Re  =  4  x  10  compared  to  the  same  projectile  with  a  frozen  liquid. 


4.  Conclusion 


An  integrated  trajectory  simulation  for  a  projectile  with  an  internal  liquid  payload  has  been 
created  and  exercised.  The  liquid  moments  are  calculated  using  linear  liquid  theory  applied  to  a 
low-viscosity  liquid  undergoing  two-mode  steady  state  coning  motion.  High  Re  boundary  layers 
approximate  the  liquid  shear  moment,  yielding  an  average  quasi-static  liquid  moment  that  is 
applied  to  a  nonlinear  6-DOF  time-dependent  trajectory  model.  These  quasi-static  averages 
yield  approximate  predictions  describing  the  motion  of  projectiles  with  low-viscosity  liquid 
payloads  filling  a  cylindrical  cavity.  The  nonlinear  6-DOF  motion  of  a  projectile  is  often  well 
approximated  as  the  sum  of  fast  and  slow  coning  motions.  Thus,  under  these  conditions  the 
linear  liquid  analysis  can  predict  when  a  projectile  exhibits  flight  instability  due  to  a  liquid 
payload  by  tracking  the  coning  frequencies  during  the  6-DOF  numerical  integration  process. 


14 


Calculations  presented  in  this  report  show  how  both  stable  and  unstable  liquid  payload 
configurations  influence  trajectory  flight  parameters.  Unstable  liquid  payload  configurations 
result  from  an  improper  combination  of  liquid  Re’s,  payload  geometry,  and  quasi-static  projectile 
coning  mode  rates.  In  such  cases,  projectile  roll  rate  and  axial  moment  decrease  rapidly  due  to 
the  quasi-static  liquid  moment. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


6-DOF 

a 

AOA 

c 

Cdd 

CFD 

I 


m 


Re 


t 

T 


six  degrees  of  freedom 

radius  of  cylinder  containing  fluid  (ft) 

angle  of  attack 

half  length  of  the  cylinder  containing  N  subcylinders  (ft) 

aero-coefficients 

computational  fluid  dynamics 

I.X  0  O' 

projectile  moment  of  inertia  I=1  0  IY  0  [  (slug  fr) 

°  0  IY 

projectile  mass  (slug) 
liquid  mass  mL  =  2  7ia2  cp  (slug) 

projectile  moment  (ft  lbf) 

projectile  angular  velocity  (s  ’1) 

displacement  of  d  relative  to  mass  center  coordinates  (ft) 

Reynolds  number,  a2p/// 
gyroscopic  stability  factor 
time  (s) 

nondimensionalized  coning  frequency 
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projectile  force  (lbf) 

nondimensionalized  growth  rate  per  cycle 
projectile  body  Euler  angles 
kinematic  viscosity  of  liquid  (ft  s  ) 

_ o 

liquid  mass  density  (slug  ft  ) 
variable  in  nonrolling  system 
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